ABSTRACT
Introduction: The B.1.1.529 (Omicron) SARS-CoV-2 variant has raised global concerns due to its high number of mutations and its rapid spread. It is of major importance to understand the impact of this variant on the acquired and induced immunity. Several preliminary studies have re-ported the impact of antibody binding and to this date, there are few studies on Omicron CD8+ T-cell immune escape. Methods: We first assessed the impact of Omicron and B.1.617.2 (Delta) variant mutations on the SARS-CoV-2 spike epitopes submitted to the Immune Epitope Database (IEDB) with positive outcome on MHC ligand or T-cell assays (n=411). From those epitopes modified by a mutation, we found the corresponding homologous epitopes in Omicron and Delta. We then ran the netMHCpan computational MHC binding prediction on the pairs of IEDB epitopes and matching homologous epitopes over top 5 MHC I alleles on some selected populations. Lastly, we applied a Fisher test to find mutations enriched for homologous epitopes with decreased predicted binding affinity. Results: We found 31 and 78 IEDB epitopes modified by Delta and Omicron mutations, respectively. The IEDB spike protein epitopes redundantly cover the protein sequence. The WT pMHC with a strong predicted binding tend to have homologous mutated pMHC with decreased binding. A similar trend is observed in Delta over all HLA genes, while in Omicron only for HLA-B and HLA-C. Finally, we obtained one and seven mutations enriched for homologous mutated pMHC with decreased MHC binding affinity in Delta and Omicron, respectively. Three of the Omicron mutations, VYY143-145del, K417N and Y505H, are replacing an aromatic or large amino acid, which are reported to be enriched in immunogenic epitopes. K417N is common with Beta variants, while Y505H and VYY143-145del are novel Omicron mutations. Conclusion: In summary, pMHC with Delta and Omicron mutations show decreased MHC binding affinity, which results in a trend specific to SARS-CoV-2 variants. Such epitopes may decrease overall presentation on different HLA alleles suggesting evasion from CD8+ T-cell responses in specific HLA alleles. However, our results show B.1.1.529 (Omicron) will not totally evade the immune system through a CD8+ immune escape mechanism. Yet, we identified mutations in B.1.1.529 (Omicron) introducing amino acids associated with increased immunogenicity.
ABSTRACT
Due to the widespread of the COVID-19 pandemic, the SARS-CoV-2 genome is evolving in diverse human populations. Several studies already reported different strains and an increase in the mutation rate. Particularly, mutations in SARS-CoV-2 spike-glycoprotein are of great interest as it mediates infection in human and recently approved mRNA vaccines are designed to induce immune responses against it. We analyzed 1,036,030 SARS-CoV-2 genome assemblies and 30,806 NGS datasets from GISAID and European Nucleotide Archive (ENA) focusing on non-synonymous mutations in the spike protein. Only around 2.5% of the samples contained the wild-type spike protein with no variation from the reference. Among the spike protein mutants, we confirmed a low mutation rate exhibiting less than 10 non-synonymous mutations in 99.6% of the analyzed sequences, but the mean and median number of spike protein mutations per sample increased over time. 5,472 distinct variants were found in total. The majority of the observed variants were recurrent, but only 21 and 14 recurrent variants were found in at least 1% of the mutant genome assemblies and NGS samples, respectively. Further, we found high-confidence subclonal variants in about 2.6% of the NGS data sets with mutant spike protein, which might indicate co-infection with various SARS-CoV-2 strains and/or intra-host evolution. Lastly, some variants might have an effect on antibody binding or T-cell recognition. These findings demonstrate the continuous importance of monitoring SARS-CoV-2 sequences for an early detection of variants that require adaptations in preventive and therapeutic strategies.